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SECTION  I 
INTRODUCTION 


In  previous  reports  (Lane,  1973;  Lane  and  Sun,l974a;  Lane  and 

Sun,  1974b,.  .he  applicability  of  the  complex  cepstrum  technique  to  the  problem 

of  separating  naturally  and  syntheticaUy  mixed  seismic  signals  into  Slmple 

waveforms  has  been  considered  for  various  cases.  Of  par.icufar  interest  in 

counter  evasion  studies  is  the  separation  of  the  signal  from  a  single  underground 

nuclear  explosion  into  direct  and  surface-reflected  phases,  and  these  studies 

showed  that  with  adequate  signal-to-noise  ratio  this  is  possible  for  signal  time 

separations  down  to  0.  1  seconds,  although  for  normally  buried  expiosions  low 

Signal-, o-noise  ratio,  limits  time  resolution  to  about  0.4  seconds.  The  then 

retie al  basis  for  the  cepstrum's  performance  was  also  presented  in  these 
papers. 

The  present  work  is  concerned  with  another  problem  in  counter- 
evasion.  that  of  the  resolution  of  multiple  underground  nuclear  explosions  into 
e.r  constituent  parts.  A  theoreticai  discussion  of  the  order  in  which  this 
s  ould  he  done  is  presented,  followed  by  die  results  of  the  cepstral  analysis  of 

a  num  er  of  seismic  events  presumed  to  be  peaceful  underground  nuclear  ex- 
plosions  (PNE's)  in  Eurasia. 
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CEPSTRAL  FILTERING  SEQUENCE 


Here  we  consider  the  problem  of  deconvolving  a  waveform  con¬ 
sisting  of  more  than  two  simple  signals.  In  such  a  case  there  are  a  number  of 
sequences  in  which  the  cepstrum  can  be  filtered.  It  will  be  shown  that  only  two 
such  sequences  are  useful  in  practice. 

Consider  the  case  of  three  signals  which  are  identical  except 
for  their  scale  factors,  taken  to  be  1,  a^  and  a^.  The  first  signal  is  repre¬ 
sented  by  s(n)  ,  and  the  signals  arrive  at  times  0,  n^,  and  n^  ,  where  the 
n.'s  are  in  units  of  the  sampling  time.  Thus  the  mixed  signal  x(n)  is 

x(n)  =  s(n)  +  a^  s(n-n^ )  +  a^n-r.^)  .  ( II - 1 ) 

Then  an  extension  of  the  development  by  Lane  and  Sun  (1974b)  shows  that  the 
complex  cepstrum  of  the  mixed  signal  is 

°°  k 

x(n)  =  s(n)  +  ^  ^  ^  ^  (-l)k+1  fijn  -  k^  +  r(n2-n2)}  (II-2) 

k=l  r=0 

where  s(n)  is  the  cepstrum  of  s(n)  ,  the  double  sum  represents  the  contribu¬ 
tion  of  the  multipath  operator,  and  the  delta  function  has  its  usual  meaning. 
Points  where  the  argument  of  the  delta  function  is  zero  are  cepstral  peaks. 

Ideally  the  effect  of  any  of  the  later  arrivals  can  be  removed  by 
zeroing  the  appropriate  cepstral  peaks.  In  practice,  however,  this  cannot  be 
done,  for  the  following  reasons.  First,  the  precise  locations  of  the  cepstral 
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peaks  are  unknown  before  filtering.  Next,  real  data  peaks  are  not  isolated 
points  but  have  some  width  and  may  therefore  interfere  with  one  another. 

Finally,  for  some  combination  of  delay  times  even  sharp  peaks  will  he  atop 
one  another.  This  will  always  happen  when,  for  example,  n^  =  2n^  .  Thus 
it  is  impossible  to  remove  the  peaks  due  to  one  later  signal  without  disturbing 
those  due  to  the  others,  even  if  accurate  knowledge  of  the  n^  is  available. 

In  this  situation  the  correct  procedure  is  to  zero  every  cepstral 
peak,  retaining  only  s(n).  This  is  done  with  the  short-pass  filter,  which  sets 
every  cepstral  point  beyond  some  time  to  zero  without  disturbing  those  at  short¬ 
er  times.  These  earlier  points  will  be  due  to  s(n)  alone  if  the  filter  point 
falls  before  the  time  of  the  first  multipathe  arrival.  From  s(n)  the  original 
signal  s(n)  in  the  sequence  can  be  recovered,  and  subtraction  of  s(n)  from 
x(n)  yields  the  remaining  waveforms.  Short-pass  filtering  can  then  be  used 
to  extract  the  first  signal  of  this  new  sequence  in  the  same  manner  as  above. 
This  procedure  can  be  repeated  a  number  of  times  but  is  ultimately  limited  to 
about  six  arrivals  by  small  but  cumulative  inaccuracies  in  the  estimation  pro¬ 
cedure. 

There  is  a  case  where  a  different  filtering  sequence  is  neces¬ 
sary.  however.  The  n-th  point  of  the  input  trace  must  be  multiplied  by  an 
exponential  weighting  factor  a°  ,  where  a  <  1.0,  in  order  that  the  multipath 
operator  be  minimum  phase  and  thus  completely  represented  in  the  cepstrum. 
The  signals  recovered  by  filtering  are  then  unweighted  by  multiplying  by  a  n. 

If  the  trace  is  long  and  a  is  small,  large  oscillations  will  appear  in  the  un¬ 
weighted  outputs  towards  their  ends  due  to  the  multiplication  of  noise  by  large 
values  of  the  factor  a~n.  These  oscillations  completely  obscure  the  motion 
of  interest,  and  can  only  be  eliminated  by  increasing  the  weighting  factor  a  , 
at  the  expense  of  making  the  multipath  operator  mixed  phase.  Useful  results 
can  still  be  obta.ned  in  the  mixed  phase  case  by  filtering  later  in  the  cepstrum, 
because  later  arrivals  may  be  represented  completely  at  large  values  of  a 


where  early  ones  are  not.  Some  distortion  will  be  introduced  by  zeroing  later 
cepstral  peaks  due  to  the  second  arrival,  so  this  sequence  is  to  be  avoided 

when'  ver  possible. 
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SECTION  III 

EXPERIMENTAL  RESULTS 


Her*'  we  present  the  results  of  cepstral  analysis  of  a  number  of 
presumed  underground  nuclear  explosions  away  from  normal  test  sites  in  Eur¬ 
asia,  assumed  to  be  for  some  peaceful  use.  All  such  events  occuring  through 
1973  for  which  good  array  data  were  available  were  examined.  Nordyke  (1974) 
presents  information  for  these  events  as  well  as  for  others  for  which  waveforms 
were  not  available.  It  is  hoped  that  these  analyses  will  prove  useful  to  others 
investigating  these  events,  and  tha*  the  results  themselves  will  help  to  estab¬ 
lish  the  validity  of  the  method. 

Fifteen  events  were  recorded  with  good  signal-to-noise  ratio  at 
LASA  during  this  time.  Three  of  them  were  far  enough  from  NORSAR  to  pro¬ 
duce  useable  teleseismic  records  there,  and  two  were  recorded  at  Mould  Bay, 
Canada.  Table  III- 1  presents  source  information  for  these  events,  and  establishes 
a  system  of  nomenclature  for  them.  Event  locations  are  shown  in  Figure  III- 1 . 

For  each  event,  the  waveform  recorded  at  LASA  was  separated 
into  its  constituent  parts,  which  were  identified  as  being  either  direct  p  phases 
or  as  surface  reflections  on  the  basis  of  the  sign  of  their  first  motion.  The 
correlation  coefficient  between  the  first  direct  and  depth  phases  was  calculated 
wherever  these  phases  were  observed,  except  for  PNE/902/69  and  PNE/908/69, 
for  which  all  analysis  was  conducted  during  a  previous  contract  period.  The 
results  are  presented  in  Table  III-2,  where  all  arrival  times  and  amplitudes 
of  later  phases  are  measured  relative  to  the  first-arriving  phase.  The  symbol 
s^(n-n.)  stands  for  a  direct  phase  delayed  by  n.  sample  points  of  0.  1  second 
duration  each.  Depth  phases  are  represented  by  s)  (n-n.).  Complicated 
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FIGURE  III- 1 


LOCATIONS  OF  PRESUMED  UNDERGROUND 
NUCLEAR  EXPLOSIONS 


later  arrivals  which  clearly  consist  of  more  than  one  phase  but  could  not  be 
resolved  with  precision  up  to  our  usual  standards  are  represented  by 
x^(n-m)  -r  by  xj  (n-m)  ,  if  the  sign  of  their  first  motion  was  positive  or  neg¬ 
ative,  respectively.  Waveforms  of  the  original  mixed  signals,  x(n),  and  those 
of  the  cepstrum  resolved  signals,  s.(n),  s!(n),  x^(n),  and  x!(n)  are  given  in 
Appendix  A. 

The  most  striking  leatures  of  this  table  are  the  large  amplitudes 
of  some  depth  phases.  The  boundary  conditions  on  stress  at  the  free  surface 
require  that  the  ratio  of  the  depth  phase  to  direct  phase  amplitude  be  -1,  and 
it  must  be  numerically  less  if  there  is  significant  spalling  or  attenuation  in 
the  near  surface  layers.  Thus  amplitude  ratios  near  2.0,  as  observed  he**-, 
are  physically  unreasonable  on  a  simple  model,  and  in  fact  have  not  been  ob¬ 
served  by  us  in  presumed  explosions  from  eastern  Kazakh. 

A  weak  decrease  in  correlation  coefficient  with  increasing  am¬ 
plitude  ratio  is  seen  in  the  data  as  plotted  in  Figure  III-?.,  but  is  probably  not 
statistically  significant.  The  average  of  the  correlation  coefficients  for  the 
presumed  PNE's  of  amplitude  ratio  near  2.0  is  0.847,  which  compares  well 
with  the  average  correlation  coefficient  of  0.848  for  13  depth  phases  from  east¬ 
ern  Kazakh  presumed  explosions.  These  measurements  indicate  that  the  de- 
cc  lpositions  with  anomalously  large  amplitude  ratio  are  nevertheless  of  good 
quality. 

Table  III—  3  shows  the  results  of  cepstral  analysis  of  the  events 
which  were  available  at  more  than  one  station,  using  the  same  notation  as  in 
Table  III-2.  Two  Kazakh  events  are  also  included  for  comparison. 

The  first  two  presumed  PNE's  show  excellent  agreement  be¬ 
tween  the  results  of  analyses  on  LASA  and  Mould  Bay  data.  This  may  be  re¬ 
lated  to  the  fact  that  the  event  azimuth  is  almost  the  same  at  Mould  Bay  as  at 
LASA.  The  other  decomposition*  are  also  compatible  with  one  another. 
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FIGURE  III -2 

CORRELATION  COEFFICIENT  VERSUS  AMPLITUDE  RATIO 


III- 6 
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Asreement  for  amplifies  and  delay  times  is  good  for  PNE/815/73  The 
ond  arrival  in  PNE/828/73  a.  NORSAR  is  a  direct  one,  a,  the  same  time  as  the 
second  direct  phase  at  LASA  and  contains  subsequent  arrivals  which  could  not 
be  resolved  due  to  their  low  level  compared  to  the  coda  from  the  first  event 
Pfre/919/73  as  recorded  at  NORSAR  also  contains  an  unresolved  later  phaee. 
whose  first  motion  is  clearly  negative,  and  which  might  be  composed  of  a 
depth  phase  plus  a  later  direct  arrival  as  in  the  LASA  analysis. 

Those  PNE  events  for  which  depth  phases  were  found  have  rea¬ 
sonable  amplitude  ratios  at  NORSAR.  when  observed  there,  of  about  half  the 
value  observed  at  LASA.  It  is  also  interesting  that  the  second  event  claimed 
for  PNE/828/73  has  a  different  amplitude  ratio  at  NORSAR  than  at  LASA. 

These  results  stregthen  the  belief  that  the  arrivals  presented 
in  Table  111-2  represent  real  seismic  waves  and  are  not  the  result  of  forcing 
the  data  to  agree  with  a  pre-assumed  model.  However,  the  amplitude  esti- 
mates  of  these  arrivals  are  quite  variable. 

Support  for  this  view  also  comes  from  the  Kazakh  events.  Here 
the  observed  delay  times  agree  within  a  sample  point,  but  the  amplitude  ratios 

disagree  by  about  30%  even  for  the  clear  high  signal-to-noise  ratio  events 
studied. 

Since  the  direct  and  reflected  phases  follow  the  same  seismic 
path  except  near  the  shot  point.  1,  1,  reasonable  to  look  there  for  an  explana¬ 
tion  of  the  observed  anomalies.  There  are  two  possibilities  -  either  the  dir¬ 
ect  phase  is  unusually  small,  or  the  surface  phase  is  unusually  large.  The 
firs,  would  arise  if  some  special  conditions  a,  the  shot  point  inhibited  the  pro¬ 
pagation  of  energy  downward.  This  might  be  the  case  since  presumably  these 
explosions  had  some  technological  purpose.  This  same  mechanism  might  also 
be  directional,  explaining  the  variation  in  amplitude  of  the  second  direct  ar¬ 
rival  of  PNE/828/73  as  well  as  the  variations  in  depth  phase  amplitudes. 
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The  second  possibility  might  occur  if  special  topography  focused 
the  surface-reflected  energy  somewhat.  The  mountainous  region  where  these 
events  occurred  might  provide  such  topography  which  might  also  focus  energy 
differently  in  different  directions.  Both  of  these  explanations  are  highly  spec¬ 
ulative,  however,  and  the  only  conclusions  drawn  here  are  that  depth  phase  to 
direct  phase  amplitude  ratios  are  quite  variable  from  event  to  event,  and  from 
recording  site  to  recording  site,  while  identification  of  phases  and  measure¬ 
ments  of  delay  times  are  probably  accurate. 

Depth  phase  delay  times  are  plotted  against  event  mb  in  Figure 
III- ^  for  the  presumed  PNE's  studied  here.  There  is  a  weak  increase  in  de¬ 
lay  time  wit  \  magnitude,  and  the  observed  times  are  not  in  conflict  with  I  'ose 
predicted  by  the  cube-root  scaling  law  shown  as  a  solid  line.  This  semi- 
empirical  rule  requires  that  contained  explosions  be  burried  at  a  depth  in  feet 
of  at  least  h  =  W  /400  ,  where  W  is  the  yield  in  kilotons  of  TNT.  An 
estimate  of  the  near-surface  p-wave  velocity  v  leads  to  a  delay  U.  le  of 
2h/v,  and  this  velocity  was  chosen  to  be  3.2  km/ sec  here.  There  is  no  rea¬ 
son  to  suppose  that  these  explosions  were  burried  near  the  depth  necessary 
for  containment,  so  divergences  of  the  data  from  the  solid  line  have  little  sig¬ 
nificance,  except  to  suggest  that  these  were  not  cratering  explosions.  The 
possible  exception  is  PNE/  1222/71,  which  fell  below  the  cube  root  scaling 

curve  at  m  =  6.  0. 
b 

Although  the  evidence  presented  here  supports  the  assumption 
that  these  events  were  multiple  explosions  it  might  possibly  be  due  to  surface 
effects  such  as  reverberations.  It  would  therefore  be  very  desirable  to  anal¬ 
yze  them,  or  others,  at  as  many  stations  as  possible.  Identical  decomposi¬ 
tions  at  three  or  more  stations  for  the  same  events  would  be  much  stronger 
evidence  of  their  nature,  especially  if  wide  azimuthal  coverage  was  obtained. 


SECTION  IV 
CONC  LUSIONS 


Fifteen  presumed  peaceful  nuclear  explosions  and  two  presumed 
Kazakh  explosions  have  been  processed  here  using  the  complex  cepstrum  tech¬ 
nique  developed  previously  by  these  authors.  Depth  phases  were  found  for  all 
but  three  of  the  PNE's,  and  could  have  been  used  as  discriminants  if  required. 
Ten  of  the  PNE's  appear  to  be  multiple  explosions,  with  delay  times  up  to  2.5 
seconds  in  the  case  of  one  possible  triple  explosion. 

The  observed  delay  times  of  the  depth  phases  are  consistent 
with  all  of  the  events  being  contained.  In  this  context  it  should  be  noted  that 
the  absence  of  a  depth  phase  does  not  lead  to  the  conclusion  that  a  crater  was 
formed,  since  the  shock  wave  responsible  for  the  depth  phase  arrives  at  the 
free  surface  while  the  latter  is  still  intact.  The  crater  is  formed  by  the  later 
expansion  of  the  cavity  gasses.  Thus  even  a  cratering  event  should  produce  a 
depth  phase. 

Correlation  between  observed  direct  and  depth  phases  was  high, 
and  for  the  cases  where  the  waveforms  could  be  observed  at  two  locations, 
agreement  between  times  of  depth  phase  arrival  was  good.  However,  anom¬ 
alously  high  amplitudes  were  observed  for  the  depth  phase  of  some  events  at 
LASA,  though  not  at  NORSAR,  Even  for  events  with  reasonable  pP  amplitudes 
at  each  array  there  was  variation  by  a  factor  of  two  among  the  ratio  between 
depth  and  direct  phases.  These  results  lead  to  the  conclusions  that: 

•  Cepstrum  analysis  can  recover  depth  phases  frc  "n  real  under¬ 

ground  nuclear  explosions  in  regions  of  complicated  geology 
in  a  large  majority  of  cases. 
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Cepstrum  analysis  is  believed  to  be  able  to  resolve  possible 
multiple  arrivals  in  the  real  teleseismic  short-period  P-waves 
of  a  presumed  explosion.  However,  the  resolved  multiple  ar¬ 
rivals  may  come  from  a  real  multiple  explosion  or  may  be  at¬ 
tributed  to  the  effect  of  surface  reverberations. 

Where  the  depth  phase  can  be  identified  its  delay  time  can  be 

determined  well  within  the  accuracy  reauired  for  discrimina¬ 
tion. 

The  cepstrum  resolved  depth  phase  amplitudes  may  be  larger 
than  the  theoretically  expected  value  by  a  factor  of  two.  This 

may  be  caused  by  the  surface  focusing  or  directioial  propaga¬ 
tion. 
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APPENDIX  A 

WAVEFORMS  FROM  THE  CEPSTRUM  ANALYSIS 

The  waveforms  of  the  original  mixed  signals  and  the  cepstrum 
resolved  signals,  which  are  represented  by  x(n),  s.(n),  s!(n),  x.(n)  and  x|(n) 

*n  Table  III-2,  are  shown  in  the  following  figures.  In  these  figures,  the  division 
along  the  time  axis  is  0.  1  second.  For  each  figure,  the  input  signal  to  the  cep¬ 
strum  analysis  is  given  first,  and  the  cepstrum  resolved  waveforms  follow. 
These  output  signals  of  the  cepstrum  analysis  arc  given  with  the  first  arrivals 
shown  as  solid  lines  and  the  second  arrivals  shown  as  dashed  lines.  For  the 
mixed  signal  which  consists  of  only  two  single  arrivals,  such  as  the  event 
PNE/  1022/71  in  Table  III-2,  the  signal  input  to  the  cepstrum  analysis  is  the 
original  mixed  signal  x(n)  and  the  output  of  the  cepstrum  analysis  is  the  final 
decomposition  shown  in  Table  III-2.  For  the  mixed  signal  which  consists  of 
more  than  two  single  arrivals,  such  as  the  event  PNE/709/72  in  Table  III-2 , 
the  input  to  the  cepstrum  analysis  can  be  either  the  original  mixed  signal  x(n) 
or  the  reduced  mixed  signal  x^(n)  (or  x|(n)).  The  final  decomposition  can  be 
obtained  by  summing  up  all  the  intermediate  outputs  of  the  cepstrum  analysis 
for  the  same  event.  As  an  example,  consider  event  PNE/709/72.  The  inter¬ 
mediate  outputs,  shown  in  Figure  A-4,  are: 

x(n)  =  s^n)  +  xj(n-5)  (A-l) 

xj(n)  =  s^n)  +  x2(n-5)  (A-2) 

x2(n)  =  s2(n)  +  s£(n-6)  .  (A-3) 


The  final  decomposition  is  obtained  by  summing  up  equations  (A-l)  through 
(A-3),  with  appropriate  time  shifts  In  equations  (A-2)  and  (A-3)  as  follows: 

x(n)  =  s  j (n)  +  xJ(n-5) 

with 

xj(n-5)  =  s  j (n- 5)  +  x2(n-10) 

and 

x2(n-10)  =  s2(n-10)  +  s^(n-16)  . 

Therefore, 

x<n)  =  Sj(n)  +  sj{n-5)  +  s2(n-10)  +  s^(n-16)  . 


A-2 


FIGURE  A- 7  a 
PNE/  1003/72 


FIGURE  A-7b 
PNE/ 1003/72 


